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A CLOSER VIEW OF THE RADIO-FIR CORRELATION: DISENTANGLING THE CONTRIBUTIONS OF 

STAR FORMATION AND AGN ACTIVITY 

I. Moric 1 2 , V. Smolcic 3,4,2 , A. Kimball 51 ', D. A. Riechers 2,7 , Z Ivezic 5 , N. Scoville 2 

ABSTRACT 

We extend the Unified Radio Catalog, a catalog of sources detected by various (NVSS, FIRST, 
WENSS, GB6) radio surveys, and SDSS, to IR wavelengths by matching it to the IRAS Point and 
Faint Source catalogs. By fitting each NVSS-selected galaxy's NUV-NIR spectral energy distribution 
(SED) with stellar population synthesis models we add to the catalog star formation rates, stellar 
masses, and attenuations. We further add information about optical emission line properties for 
NVSS-selected galaxies with available SDSS spectroscopy. Using an NVSS 20 cm (Fi.4GHz > 2.5 mJy) 
selected sample, matched to the SDSS spectroscopic ("main" galaxy and quasar) catalogs and IRAS 
data (0.04 < z < 0.2) we perform an in depth analysis of the radio-FIR correlation for various 
types of galaxies, separated into i) quasars, ii) star forming, iii) composite, iv) Seyfert, v) LINER 
and vi) absorption line galaxies using the standard optical spectroscopic diagnostic tools. We utilize 
SED-based star formation rates to independently quantify the source of radio and FIR emission 
in our galaxies. Our results show that Seyfert galaxies have FIR/radio ratios lower than, but still 
within the scatter of, the canonical value due to an additional (likely AGN) contribution to their radio 
continuum emission. Furthermore, IR-detected absorption and LINER galaxies are on average strongly 
dominated by AGN activity in both their FIR and radio emission; however their average FIR/radio 
ratio is consistent with that expected for star forming galaxies. In summary, we find that most AGN- 
containing galaxies in our NVSS-IRAS-SDSS sample have FIR/radio flux ratios indistinguishable from 
those of the star-forming galaxies that define the radio-FIR correlation. Thus, attempts to separate 
AGNs from star-forming galaxies by their FIR/radio flux ratios alone can separate only a small fraction 
of the AGNs, such as the radio-loud quasars. 

Subject headings: galaxies: fundamental parameters - galaxies: active, evolution - cosmology: obser- 
vations - radio continuum: galaxies 



1. INTRODUCTION 



The radio-FIR correlation is one of the tightest cor- 
relations in observat i onal astrophysics fe.g. iHelou et alJ 



1985; Condon 199 21: iMauch fc Sadler 1 120 07: Yu n et al l 
200l[|Bell 2003; Sa rgent et al. ll2010l:lKovacs et al.l l20'oa 
Murp hy et all l2009UAppleton et al. 1 120041 ) . The corre- 



spondence between the radiation in the (far-)infrared and 
that in the radio spans over nearly five orders of magni- 
tude in various types of galaxies, ranging from dwarfs to 
ULIRGs. Given that the two observational windows, IR 
and radio, trace independent and different intrinsic phys- 
ical mechanisms in galaxies - thermal vs. synchrotron 
radiation - the existence of such a tight correspondence 
is remarkable. It is generally believed that recent star 
formation in galaxies is the process that relates IR and 

Based on observations with the National Radio Astronomy 
Observatory which is a facility of the National Science Founda- 
tion operated under cooperative agreement by Associated Uni- 
versities, Inc. 

1 University of Zagreb, Physics Department, Bijenifika cesta 
32, 10002 Zagreb, Croatia 

2 California Institute of Technology, MC 249-17, 1200 East 
California Boulevard, Pasadena, CA 91125 

3 ESO ALMA COFUND Fellow, European Southern Ob- 
servatory, Karl-Schwarzschild-Strasse 2, 85748 Garching b. 
Muenchen, Germany 

4 Argelander Institut for Astronomy, Auf dem Hiigel 71, Bonn, 
53121, Germany 

5 Department of Astronomy, University of Washington, Box 
351580, Seattle, WA 98195-1580 

6 National Radio Astronomy Observatory, 520 Edgemont 
Road, Charlottesville, VA 22903 

7 Hubble Fellow 



radio emission. 

The radio-FIR correlation has been extensively 
studied in the past both in th e low dHelou et al.1 1 1985 ; 
Condon I IT992I: IGarrett I [20021: IMauch fe Sadler"! 12007 



Yun et alj 120011: iBell 1 120031) and high redshift universe 
dSargent et al. 1 120101: iKovacs et al.l l2006t iSaiina et al. I 
1 20081: IMurphv et al I 120091: lAppleton et al. I 12004: 
Vlahakis et al. 1 120081: llbar et all l2008t iChapman et al I 
I2005T ). It has been shown that, out to redshifts of 
z ~ 3 — 4, the FIR/radio ratios of various types of 
galaxies are essentially the same as those in the local 
universe. At higher redshifts radio quiet QSOs have 
been demonstrated to have FIR/radio ratios consistent 
with the local value, while the FIR/radio ratios of 
z > 4 SMGs are found to be lower by a few factors. 
This is somewhat contrary to expectations, as the the 
FIR/radio ratio is expected to be rising with redshift 
(especially at z > 3) due to the increase of the cosmic 
microwave background (CMB) energy density (C/cmb) 
with redshift, Ucmb oc (1 + z) 4 , which surpresses the 
non-thermal component of a galaxy's radio continuum 
via inverse- Compton (IC) scattering (see Murphy 2009 
for details). An explanation for this discrepancy can be 
provided by additional processes that add to a galaxy's 
radio continuum, such as increased magnetic field 
strengths or AGN contribution, that may compensate 
for the radio continuum emission losses due to IC 
scattering. 

The AGN contribution to the radio-FIR correlation has 
been studied in the past to some extent. Typically a low 
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FIR/radio ratio, significantly offsetting a galaxy from the 
correlation is thought to indicate a radio-loud AGN (e.g . 
lYun et al.1 12001 ICondon. Cotton fc Broderickl 12002) . 
However, recent studies have shown that optically- 
selected AGN often follow the correlation, albeit with 
a slightly lower FIR/rad io ratio. For exam ple, based on 
SDSS-NVSS-IRAS data. lObric etail (f2006h have demon- 
strated a tight correlation between radio and 60 fim. 
fluxes for low-luminosity AGN (predominantly Seyferts 
and LINERs), which varies by only ~ 20% relative 
to that of star fo rming galaxies. Utiliz ing 6dFGS- 
NVSS-IRAS data, IMauch fc Sadlerl l[2007ft inferred a 
lower average FIR/radio ratio for AGN-bearing galax- 
ies (Seyferts, LINERs, and quasars), but still within the 
scatter of the correlation for star forming galaxies. Fur- 
thermore, studies of the correlation at higher redshifts 
have yielded a handful of interesting objects for which it 
has clearly been shown that a significant AGN contribu- 
tion to IR and/or radio exists, yet their FIR/radio ratio 
is consistent with the canonical value for star forming 
galaxies (Riechers et al. 2009, Murphy et al. 2009). 

In order to understand in more detail the contribu- 
tion of AGN activity to the radio-FIR correlation, we 
perform an in-depth study of the radio-FIR correlation, 
with a large sample, as a function of galaxy type, and 
comparison with star formation rates for those individ- 
ual samples. The various types of star forming and 
AGN-bearing galaxie s have b een drawn from the NVS S 
(ICondon et al. 1 119981) IRAS (|Neugebauer et al. Ill98l . 
and SPSS (jYork et al]l20"00h sky surveys. In Sec. [2] we 
present the data used in this paper. We present the corre- 
lation for various types of galaxies in Sec.[3J In Sec. 2] we 
link the FIR and radio emission from galaxies in our sam- 
ple to independently derived star formation rates, and in 
Sec. [5]and Sec. [6] we discuss and summarize our results, 
respectively. We adopt H = 70, Qm = 0.3, fl\ = 0.7, 
and define the radio synchrotron spectrum as F v oc v~ a . 
assuming a — 0.7. Throughout the text we will often 
use the term 'quasar' referring to both quasi-stellar ra- 
dio sources and quasi-stellar objects. 

2. DATA AND GALAXY SAMPLES: EXPANDING THE 
UNIFIED RADIO CATALOG 

2.1. Unified Radio Catalog 

Kimball ~fe Ivezid (|2008l ) have constructed 
log of r adio sources d etecte d by the GB6 (6 cm) 
FIRS T (|Becker et al. I 11993) . NVSS (ICondon et al. I 
Il998t 20 cm), and WENSS (92 cm) radio surveys, a s 
well as the SDSS (DR6) optical survey (|York et al.ll2000T ). 
This "Unified Radio Catalog" has been generated in such 
a way that it allows a broad range of 20 cm based sample 
selections and source analysis (see Kimball & lvezidl2008l 
for details). The 2.7 million entries are comprised of the 
closest three FIRST to NVSS matches (within 30") and 
vice-versa, as well as unmatched sources from each sur- 
vey. All entries have been supplemented by data from the 
other radio and optical surveys, where available. Here 
we select from the Unified Radio Catalog (version 1.1) all 
20 cm sources that have been detected by the NVSS radio 
surve y (using matchflag_nvss = —1 and matchflagjirst < 
1; see Kimball fc Ivezicll2008l for details). This selection 
yields a radio flux limited (-Fi.4GHz ^ 2.5 mJy) sample 
that contains 1,814,748 galaxies. In the following section 
we expand this catalog to IR wavelengths, and augment 



it with additional (spectroscopic and SED-based) infor- 
mation. 

2.2. Expanding the Unified Radio Catalog 
2.2.1. IRAS 

For the purpose of this paper, we have expanded 
the Unified Radio Catalog to IR wavelengths by cross- 
correlating it with the IRAS point-source and faint- 
source catalogs (hereafter PSC and FSC, resp.). The 
IRAS PSC contains 245,889 confirmed point sources 

detected at 1 2, 25, 60 and 100 /Ltm, respectively 

iStrauss et al. I (|1990f ). The completeness of the catalog 
at these wavelengths reaches down to 0.4, 0.5, 0.6 and 
1.0 Jy, respectively. The FSC was tuned to fainter levels 
based on the same IR data by point-source filtering the 
individual detector data streams and then coadding those 
using a trimmed-average algorithm (see iMoshir et al. I 
11992ft . The reliability of the FSC is slightly lower than 
that of the PSC (> 94% compared to 99.997%); however 
its sensitivity is higher by a factor of ~ 2.5. The FSC 
contains 173,044 point sources with flux densities typi- 
cally greater than 0.2 Jy at 12, 25 and 60 ^im and 0.5 Jy 
at 100 lira. 

We used match ing radius of 30" , as optimized by 
IQbric et~a l. (2006), in cross-corelating the Unified Ra- 
dio Catalog with the IR IRAS data. In Fig. [T] we show 
the distribution of the distances between the IR and ra- 
dio detections. The cumulative distribution displayed in 
Fig. Q] shows that ~ 70% of the positional matches are 
within an angular distance of 15". 

Our NVSS-selected radio sample contains 18,313 
galaxies with high quality IR photometrj0 (see 
table:class ). As the FSC and PSC have been gener- 
ated based on the same data, most of the PSC sources 
are included in the FSC. In our entire NVSS-IRAS sam- 
ple, 26% of the sources have a PSC detection but are 
not included in the FSC. This fraction, however, reduces 
to only 3% after an optical (SDSS) cross-match is per- 
formed. 

The 60 and 100 /j,m magnitudes reported in the PSC 
and FSC are in agreement for the union of the two IR 
samples. The biweighted mean of the flux difference (for 
a subsample with SDSS detections) is 0.02 and 0.03 Jy 
at 60 and 100 /im, respectively. The root mean scatter 
of the 60 /im flux difference distribution is 0.06 Jy, while 
that of the 100 /im distribution is significantly larger, i.e. 
0.16 Jy. Therefore, in order to access the highest qual- 
ity IR photometry, hereafter we use the values reported 
in either the FSC or PSC catalog corresponding to the 
higher photometric quality flag quoted in the catalogs. 
The distribution of the 60 /im and 20 cm flux densities 
are shown in Fig. [2] 

2.2.2. SDSS quasar and main galaxy sample catalogs 

We have further matched the NVSS-selected sample 
from the Unified Radio Catalog with data drawn from i ) 
the SDSS DR5 quasar sample ([Schneider et al. 1 120071) . 
and ii) the DR4 "main" spectroscopic sample for which 
derivations of emission line fluxes from the SDSS spectra 

9 We take the IRAS quality indicator, reported in the FSC and 
PSC, to be > 2 at 60 and 100 fim (the wavelength bands utilized 
here) . 
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are available (see Smolcic 2009 and references therein; 
note that the DR5 quasar and DR4 main galaxy cata- 
logs were the most up-to-date versions available at the 
time). The latter was complemented with stellar masses, 
star formation rates, dust attenuations, ages, metallici- 
ties, and a variety of other parameters based on spectral 
energy distribution ( SEP) fitting o f the SPSS (ugriz) 
photometry using the lBruzual et all (j2003l) stellar popu- 
lation synthesis models. The SEP fitting was performed 
as described in detail in Smolcic et al. (2008). 

During the inspection of the validity of the final cat- 
alog, we have found about 1% of objects with differ- 
ent spectroscopic redshifts in various SDSS data releases 
(Az > 5 x 10 -4 ). We have excluded those from the 
sample. Furthermore, a small number (~ 0.2%) of dupli- 
cate objects was present in both the SDSS "main" galaxy 
sample and the SDSS Quasar Catalog. Visually inspect- 
ing their spectra yielded that most of these objects are 
better matched to the properties of the "main" galaxy 
sample (as no power-law continuum nor broad emission- 
lines were present in the spectrum), and we have ex- 
cluded these from our quasar sample. A summary of the 
various radio-IR-optical samples is given in table:class , 
and in Fig. [3] and Fig. @]we show the radio (20 cm), op- 
tical (r band), and far-IR luminosities as a function of 
redshift for the final NVSS-SDSS and NVSS-SDSS-IRAS 
samples (see eqs. [3] and 2]). Note that the shallow IRAS 
sensitivity (compared to the NVSS and SDSS data) sig- 
nificantly reduces the number of objects, and biases the 
sample toward lower redshifts. 

2.3. Radio- optical- IR samples 
2.3.1. Star forming and AGN galaxy subsamples 

We have used the optical spectroscopic informa- 
tion added to the NVSS selected sample to spectro- 
scopically separate the galaxies present in the SDSS 
(DR4) "main" galaxy sample as absorption line, AGN 
(LINER/Seyfert), star forming, or composite galaxies. 

We define emission-line galaxies as those where the rel- 
evant emission lines (Ha, H/3, O[III,A5007], N[II,A6584], 
S[II,AA6717,6731]) have been detected at S/N > 3, and 
consider all galaxies with S/N < 3 in any of these 
lines as absorption line systems (see e.g . iBest et al"1l2005t 
iKewlev et al.ll2006HSmolcic et al.ll2009l ). As strong emis- 
sion lines are not present in the spectra of the latter, 
yet they are luminous at 20 cm, they can be consid- 
ered to be (low excitati on) AGN (see e.g. IBest et all 
2005; ISmolcic et aD[2008l for a more detailed discussion) . 
Further, as illustrated in Fig. El using standard optical 
spectroscopic diagnostics ([Baldwin. Phillips fc Terlevichl 
U98lHKauffmann et al.ll2003bl; IKewlev et al.ll200lL I2006D 
we sort the emission-line galaxies into 1) star forming, 
2) composite, 3) Seyfert, and 4) LINER galaxies. The 
last two classes have been selected 'unambiguously' by 
requiring combined criteria using three emission line flux 
ratios (see middle and right panels in Fig. [3]). A sum- 
mary of the number of objects in each class is given 
in tablexlass . It is noteworthy that the IR detec- 
tion fraction is a strong function of spectral class. It is 
the lowest for absorption line (0.6%) and LINER (6.5%) 
galaxies, intermediate for Seyferts (22%) and the highest 
for composite (40%) and star forming (46%) galaxies. 
These results suggest lower amounts of dust (and gas; 
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ISolomon fe Vanden Boutl [2005) in the former or alter- 
natively dominantly very cold dust that peaks at longer 
wavelenghts. 

The redshift distribution of the various galaxy types 
with 20 cm NVSS and NVSS-IRAS detections is shown 
in the two top panels in Fig. [5] Note that the redshift 
distribution of 20 cm detected absorption line galaxies 
is biased toward higher-redshifts, compared to all other 
galaxy types (see top panel in Fig. [B]). However, this 
is not the case when an IRAS IR detection is required, 
as illustrated in the middle panel in Fig. |6] The IR de- 
tection fraction of the different galaxy classes is shown 
as a function of redshift in the bottom panel in Fig. [51 
Except for the overall trend that absorption and LINER 
galaxies are detected less efficiently in the IR, there is no 
substantial difference between the detection fractions as 
a function of redshift for different types of spectroscopi- 
cally selected galaxies. 

Hereafter, we apply redshift range limits of 0.04 < z 
< 0. 3 to our sample. The lower redshift limit is adopted 
from IKewlev et all (|2005l ). Kewley et al. explored effects 
of fixed-size aperture of the SDSS spectroscopic fibers 
on the spectral characteristics such as metallicity, star 
formation rate, and reddening. They concluded that a 
minimum aperture size covering ss 20% of spectral light 
was required to properly approximate global values. The 
SDSS fiber aperture of 3" diameter collects such a frac- 
tion of light for galaxies of average size, type, and lumi- 
nosity at z > 0.04. The upper redshift limit of z = 0.3 
is equivalent to that of the SDSS "main" spectroscopic 
sample (note however that the majority of IR-detected 
galaxies are at z < 0.2, see Fig. [J}. It is worth noting 
that, because of lower spectral signal-to-noise for fixed- 
luminosity galaxies at greater distances, galaxies with 
weak emission lines, such as LINERs, ca n get confused 
with absorption line galaxies at z > 0.1 (|Kewlev et al.l 
2005). However, as LINER and absorption galaxies have 
similar physical properties, (e.g. Smolcic 2009) we sim- 
ply combine these two types of galaxies, and treat them 
hereafter as a single class. 

2.3.2. Quasar subsample 

Matching the SDSS DR5 quasar catalog to the Unified 
radio catalog resulted in 4,490 matches (see tablexlass 
). The redshift range of our radio luminous quasars 
is 0.09 — 5.12, with a median at z = 1.36. Requiring 
IRAS detections biases the sample toward low redshifts 
(0.12 < z < 1.15), with a median redshift of 0.18, and 
selects only ~ 0.5% of the radio detected quasars. The 
radio (> 10 23 W Hz" 1 ) and FIR (> 2 x 10 11 L ) lumi- 
nosities (see eqs. [3] and [4]) of our quasars are systemati- 
cally higher than those of the SDSS "main" spectroscopic 
sample galaxies in our radio-optical-IR sample. 

3. QUANTIFYING THE RADIO-FIR CORRELATION FOR 
VARIOUS SOURCE TYPES 

3.1. Parameterizing the radio-FIR correlation 

The radio-FIR correlati on is usually quan tified by its 
slope via the q parameter ()Helou et al.lll985T) , defined as 
the logarithmic ratio of the far-infrared flux to radio flux 
density: 

. f f FIR /(3.75xl0 12 Hz) ) 
9 = log<^ p > (1) 
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where -Fi.4GHz is the 1.4 GHz radio flux density in units 
of Wm~ 2 Hz _1 and Jfir. is the far-infra r ed flu x in units of 
Wm~ 2 . Following [Sanders fe Mirabell ([1996H we define 
the latter as: 

F Fm = 1.26 x 1(T 14 (2.585 60 ^m + Sioo^m) (2) 

where Seo/im and >Sioo/zm are observed flux densities 
at 60 and 100/im (in Jy), respectively. 
We compute the far-infrared luminosity as: 

L FIR = 4TrD 2 L CF FIR [LQ] (3) 

where Dl is the luminosity distance (in units of m) 
and C is a scale factor used to correct for the extrapo- 
lated flux longward of the IRAS 100 (tm filter. We use 
C = 1.6 (see Tab. 1 in lSanders fe Mirabell [1991 . Note 
that this expression can also be utilized to compute the 
FIR luminosities for our IR-detected quasars, given their 
relatively low redshifts. 

The radio luminosity density is computed as: 

M.4GHz - ^ + v 1 _ g i'l.4GHz \V 

where z is the redshift of the source, -Fi.4GHz is its in- 
tegrated flux density, and a the radio spectral index (as- 
suming F v oc v~ a ). To compute the radio luminosities, 
we assumed a spectral index of a = 0.7. 

3.2. Radio-FIR correlation for all sources 

The radio-FIR correlation for the NVSS-SDSS-IRAS 
sample is summarized in Fig. [7J The radio and FIR 
flux densities (top left panel) and luminosities (top right 
panel) clearly show a tight correlation that holds over 
many orders of magnitude. In the middle panels we 
show the q parameter, that characterizes the slope of 
the radio-FIR correlation (see eq. H}, as a function of 
FIR and radio luminosities. The average q is constant as 
a function of FIR luminosity (middle left panel), and it 
is decreasing with increasing radio power (middle right 
panel; see also below). In the bottom panels of Fig. [7] 
we show the q parameter as a function of redshift, as 
well as its distribution for all our NVSS-SDSS-IRAS 
sources (galaxies and quasars). We find that the av- 
erage (biweighted mean) q value for the entire NVSS- 
SDSS-IRAS sample is q=2.273±0.008, with a root-mean- 
square scatter of a = 0.18. This is in very good agree- 
ment with previous findings ([Condon I [1992 : lYun et all 
20fUL ICondon. Cotton fc Broderickl 120021 : IBell I 12001 
Mauch fe Sadler 1 120071 ). and will be discussed in more 
detail in Sec. 

The quasars in our sample comprise the high luminos- 
ity end at both IR and radio wavelengths (they are also 
located at higher redshifts, compared to the IR and radio 
detected "main" galaxy sample). It is also obvious that 
there is a larger fraction of quasars that do not lie on the 
radio-IR correlation, compared to that for the "main" 
sample galaxies. 

3.3. Radio-FIR correlation for different types of galaxies 

In Fig. [5] we present the radio-FIR correlation for 
the SDSS "main" galaxy sample subdivided into differ- 
ent, spectroscopically selected galaxy types (absorption, 



LINER, Seyfert, composite and star forming galaxies; see 
Sec. 12.3.11 and Fig. [5] for details on the selection). The 
decrease in q with increasing radio power for all types of 
galaxies (middle right panel) is consistent with various 
other observations (e.g. Smolcic et al. 2008; Kartaltepe 
et al. 2010; Sargent et al. 2010, Ivison et al. 2010). Note 
that, given the definition of the FIR/radio ratio, for any 
sample in which q does not vary with FIR luminosity, and 
has a non-zero dispersion, it is expected to decrease with 
increasing radio luminosity (see e.g. Condon 1984). To 
test whether the magnitude of the decrease is as expected 
from statistics or e.g. higher due to an additional effect 
(such as AGN contribution) we have computed the radio 
luminosity for each source based on its observed FIR lu- 
minosity and a FIR/radio ratio drawn from a Gaussian 
distribution with a dispersion of 0.18 and a mean of 2.27. 
We find that the decrease of q with radio luminosity in 
the observed data is consitent with that in the simulated 
data, thus not requiring additional effects (such as in- 
creasing AGN contribution with increasing radio power) 
to explain this trend (at least in the radio luminosity 
range probed here). 

A quantitative analysis of the radio-FIR correlation for 
different galaxy types is presented in Fig [9] The spectro- 
scopic selection of pure star forming galaxies allows us 
to quantify the radio-IR correlation in a rather unbiased 
manner. For our star forming galaxies we find an average 
q value of < q > = 2.27 ± 0.05, with a small root-mean- 
square scatter of a = 0.13. It is interesting to note that as 
the AGN contribution rises in galaxies (as inferred based 
on optical spectroscopic diagnostics) the scatter in q in- 
creases by ~ 50% to ~ 150%. Interestingly, the scatter 
is the highest for Seyfert types of galaxies, for which we 
also find the lowest average q value, < q > = 2.14 ± 0.05. 
These differences will further be discussed in Sec. [5] 

3.4. Radio-FIR correlation for quasars 

In Fig.[TU]we quantify the radio-FIR correlation for the 
21 IR-detected quasars in our sample. The distribution 
of the FIR/radio ratio cannot be well fit with a Gaus- 
sian distribution. The median q value of the sample is 
2.04, comparable to the average q value we have found 
for Seyfert galaxies (2.14), and lower than that for star 
forming galaxies (2.27; see Fig. It is worth noting 
that the higher redshift quasars (0.2 < z < 0.4) appear 
to be biased toward more radio-loud AGN. 

4. AN INDEPENDENT VIEW OF THE RADIO-FIR 
CORRELATION: A LINK TO STAR FORMATION 

It is generally taken that recent star formation 
drives both the radio and FIR emiss ion of galaxies 
that lie on the radio- FIR correlation ((Condon 1 119921 : 
iMauch fe Sadlerl [20071 Therefore, a correlation is ex- 
pected to be present between the SFRs and radio/FIR lu- 
minosities obtained from the fluxes of galaxies dominated 
by recent star formation. To shed light on the source of 
radio/FIR emission in our galaxies, in this section we in- 
vestigate the correlation between their radio/FIR lumi- 
nosities and star formation rates, independently deter- 
mined based on fitting stellar population synthesis mod- 
els to the NUV-NIR SED). 

We have derived a star fo rmation rate for every galaxy 
in our sample by fitting the iBruzual et all ([20030 library 
of stellar population synthesis models to the SDSS ugriz 
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photometry (see Sec. H}. In Fig. QT] we show the ra- 
dio and FIR luminosities of our sources as a function 
of our SED-based SFRs. As expected, a correlation is 
discernible between these two quantities. This is espe- 
cially emphasized for star formation dominated galaxies 
(i.e. star forming and composite galaxies), the distribu- 
tion of which agrees well with the co mmonly used ra- 
dio /IR luminosi ty - SFR calibrations ()Kennicutt 1 119981 : 
lYun et all 120011 ). Note that this is quite remarkable as 
the SFRs have been derived completely independently 
from the FIR or radio emissions in the galaxies. 

From Fig.[TT]it is obvious that a large fraction of galax- 
ies with significant AGN contribution (Seyfert, LINER 
and absorption galaxies) has an obvious excess of radio 
power and FIR luminosity compared to that expected 
from the galaxy's SFR. The most obvious example of 
this are the LINER and absorption galaxies from both 
the NVSS-SDSS and NVSS-SDSS-IRAS samples. 

To investigate whether star formation is the under- 
lying source of radio/FIR emission in our galaxies, we 
further quantify the difference between the radio/FIR 
emission and that expected from star formation. We 
thus define an "excess" in radio/FIR emission relative 
to that expected from star formation ( A log L\ ag\o, and 
AlogLpiR, resp.) as: 



A log L = log L. 



data 



log L c 



(5) 



where logLdata is the logarithm of the 1.4 GHz or FIR 
luminosity derived based on NVSS or IRAS data (see 
Sec. [5]), and logL cxp . is the luminosity (either at 1.4 GHz 
or FIR) expected based on the SED-derived SFR and 
the standard radio and FIR luminosity to SFR calibra- 
tions. To convert SF R to radio lumin osity we use the 
calibration defined in lYun et all (|2001ft : SFR [M Q /yr] = 
5.9xlCr 22 Li.4GHz [W/Hz]. To convert SFR to FIR lumi- 
nosity we use the standard conversion defined by Kenni- 
cutt (1998): SFR [M Q /yr] = 4.5 x 1(T 37 L Fm [W]. Prior 
to applying these conversions, derived using a Salpeter 
IMF, we have scaled our SED-based SFRs by -0.2 dex to 
convert from a Chabrier to a Salpeter IMF (we have addi- 
tionally included a scaling factor of ~ 0.4 dex to account 
for the star formation histories used in our models; see 
Smolcic et al. 2008 and Walcher et al. 2008 for details). 

Fig. QJ shows the FIR (AlogLpm) versus 1.4 GHz 
(A log Li.4GHz) luminosity excess for different types of 
galaxies in the NVSS-SDSS-IRAS sample. As expected, 
star forming galaxies follow a normal distribution in 
both the FIR and radio luminosity excess, with a mean 
AlogL value of about zero (< AlogLpm >= 0.001, 
< AlogLi.4GHz >= 0.06). The root-mean-square scat- 
ter is 0.35, and 0.32 for the FIR and radio distributions, 
respectively. Assuming the validity of the SFR to ra- 
dio/FIR luminosity calibrations, such a (normal) distri- 
bution is expected if FIR and radio emissions arise from 
star formation processes in the galaxies. From Fig. [T^] 
it is apparent that, as the AGN contribution (defined 
via optical spectroscopic emission line properties) rises 
in galaxies, the distribution of both the FIR and radio 
luminosity excess becomes highly skewed towards higher 
A log L values. Although for a fraction of optically se- 
lected AGN it is possible that the AGN contribution to 
the radio/FIR may be weak (AlogL ~ 0), and they 
may be overwhelmed by star formation (which results 



in the canonical FIR/radio ratio), the significant skew- 
ness of the A log L distributions suggests an additional 
source of radio and FIR emission in AGN bearing galax- 
ies (at least for AlogL > 0). Even more interesting is 
that galaxies with large luminosity excess in both FIR 
and radio emission (AlogL > 0.9) predominantly have 
FIR/radio ratios consistent with the mean q value for star 
forming galaxies (see Fig. 12). This will be discussed in 
more detail in Sec. [5] 

5. DISCUSSION 



5.1. Comparison with previous results 



of the r a dio-FI R correlatio n 
iCondon I jl992l: lYun et all 
Condon. Cotton fc Broderick I [2002; O bric et all 



Extensive studie s 
Helou et al.1 H98l 



2001 



2006 



Mauch fc Sadler 1120071 ) have led to an average FIR/radio 
ratio in the local (z < 0.3) universe of q ~ 2.3, and lower 
for A GN-bearing galaxies (see Tab. 2 in ISargent et al. I 
2010 for a summary). For example, using the IRAS 2 Jy 
galaxy sample (i^o/^m > 2 Jy; 1809 sources with op- 
tical counterparts and w ell determi n ed re dshifts) com- 
bined with NVSS data, lYun et all (|2001[ ) have found 
< q >= 2.34 ±0.01. A lower average q value is generally 
inferred when using r adio selected samples , and reaching 
fainter in the IR ( see ISargent et al. 1 120101 for a detailed 
discussion of selection effects). Combining NVSS data 
with the optical Uppsal a Galaxy Catalog (UGC) and 
the IRAS FSC and PSC, ICondon. Cotton fc Brodenckl 
(2002|) have found < q >= 2.3 and rms width a = 0.18. 
Furthermore, matching N VSS and 6dFGS s u rvey data 
only with the IRAS FSC, IMauch fc Sadlerl (f2007h in- 
ferred a mean q value of 2.28 with a root-mean-square 
scatter of 0.22 for their entire sample. For a subset of 
radio-loud AGN (that would correspond to our Seyfert, 
LINER, absorption, and quasar classes combined) they 
found an even lower average value, < q >= 2.0, and a 
sigificantly higher scatter in the FIR/radio ratio (a — 
0.5). 

In Fig. [9] we have presented the distribution of q for 
various types of our spectroscopically selected NVSS- 
SDSS-IRAS (PSC+FSC) galaxies. Our results yield that 
the dispersion is the tightest for star forming galaxies 
(a = 0.13), and rises by a factor of 1.5, 2.5, and 2.2 for 
composite, Seyfert and absorption/LINER galaxies, re- 
spectively. We find that the average FIR/radio ratio for 
all objects in our radio-optical-IR sample is 2.27 ± 0.01 
with a dispersion of .2. This is in very good agreement 
with the results from IMauch fc Sadler"! (|2007f). Further- 
more, if we limit the 60 /im fluxes of our full sample to 
> 2 Jy we obtain a n average value o f 2.34, consistent 
with that inferred by lYun et all (|2001l ) . 

Our results yield a lower FIR/radio ratio (< q > = 
2.14±0.05) for Seyfert galaxies, and a significantly higher 
root- mean-square-scatter (a — 0.3), compared to that 
found for SF galaxies. It is interesting that the mean 
q value for our IR-detected LINER and absorption line 
galaxies is comparable to that for star forming galaxies. 
However, the spread in q for the former is significantly 
larger than for the latter (0.28 compared to 0.13, respec- 
tively). The average FIR/radio ratio for the 21 quasars 
in our sample is q = 2.04, comparable to that inferred for 
Seyferts and lower than that for star forming galaxies. If 
we combine our AGN-bearing galaxies (quasars, Seyferts, 
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LINERs, absorption galaxi es) into one class i n ord er to 
match the AGN sample of lMauch fe Sadler I (|2007t ). we 
infer an average q of 2.16±0.03 (with a root-mean-square 
scatter of a = 0.24). This is in relatively good agreement 
with their results. In the next sections we will discuss the 
variation of q with radio luminosity and the AGN con- 
tribution to the radio-FIR correlation. 

5.2. AGN contribution to the radio-FIR correlation 

A low q value is often used to discriminate be- 
twe en star forming galaxies and AG N. For exam- 
ple iCondon. Cotton fe B rodcrick (2002) have classified 
radio-loud AGN as those having q < 1.8. Assuming that 
the FIR emission arises solely from star formation, this 
criterion selects galaxies with more than 3 times the ra- 
dio emission from galaxies on the FIR-radio correlation. 
lYun et aTl ()2001h have used q = 1.64 as a star forma- 
tion/AGN separator, identifying galaxies that emit in 
radio more than 5 times than predicted by the correla- 
tion. It is important to point out that these discriminat- 
ing values are tuned to select only the most radio-loud 
AGN. Having i) separated our NVSS-SDSS-IRAS sample 
into various classes of AGN, and ii) independently esti- 
mated SFRs in their host galaxies, we can now analyze 
the physical source of FIR and radio emission in galaxies 
both following and offset from the radio-FIR correlation. 

Assuming that the additional source of FIR and ra- 
dio emission (relative to that expected from star forma- 
tion) observed in our composite, Seyfert, absorption and 
LINER galaxies (see Fig. [T^] ) arises from the central su- 
permassive black hole, the distribution of our luminosity 
excess, AlogL defined in eq. [5j allows us to constrain 
the average contribution of star formation and AGN ac- 
tivity to the total power output for a given galaxy pop- 
ulation. Taking that star formation and AGN activity 
are the two dominant FIR/radio emission generators, 
i.e. L to t = ^sf + -^agn, the average fractional contri- 
butions of these two sources (< /sf >, and < /agn >) 
to the total power output can then be computed as 
< /sf >= 10- <A1 °s i> , and < /agn >= 1- < /sf >, 
where < A log L > denotes the average (median) of the 
AlogL distribution (see Fig. IT2"1). 

The median AlogL values, and the fractional 
star formation/ AGN contributions are summarized in 
table:fracs . As expected, for star forming galaxies we 
infer that the average contribution to FIR and radio 
emission due to star formation is 100%. We find that 
composite objects are dominated by star formation at 
the ~ 80 — 90% level. Further, the FIR emission from 
Seyfert galaxies arises predominantly from star formation 
(~ 75%), while the AGN contribution to radio luminos- 
ity in Seyfert galaxies is about a factor of two higher in 
the radio than in the FIR (see table:fracs ). The latter 
explains the lower average q value (compared to the nom- 
inal ^yjahjBjJjoj^eyfert ^ as in 
e.g. lObric et all (l2006T ). [Mauch fc Sadler I (f2007h . Lastly, 
based on the above calculation IR detected absorption 
and LINER galaxies are on average strongly dominated 
by AGN activity (~ 90%) in both their FIR and radio 
emission although their average FIR/radio ratio is con- 
sistent with that expected for star forming galaxies (see 
Fig. IS). 

One of the main results of this work is that, for the 



large majority of galaxies with radio and/or IR emis- 
sion excess, we infer q values consistent with the aver- 
age FIR/radio ratio found for star forming galaxies (see 
Fig. fT2)) . Thus, although a significant AGN contribution 
is likely present in these galaxies (adding both to the 
FIR and radio emission), they would not be identified 
with a simple low-q discriminator value, as is commonly 
used to select radio-loud AGN. Our results indicate that 
the FIR/radio ratio is not particularly sensitive to AGN 
contribution and that the radio-FIR correlation is a poor 
discriminant of AGN activity, except for the most power- 
ful AGN. This is consistent with observations of several 
AGN-bearing galaxies in the high-redshift universe. 

Based on an SED analysis, Riechers et al. (2009) find 
that both the radio and FIR luminosity in the z=3.9 
quasar APM08279+5255 are dominated by the central 
AGN, but that it has a q value consistent with the lo- 
cal radio-FIR correlation. Furthermore, Murphy et al. 
(2009) have analyzed Spitzer-IRS spectra of a sample of 
22 0.6 < z < 2.6 galaxies, composed of submillimeter 
galaxies, as well as X-ray and optically selected AGN in 
GOODS-N. Making use of their IRS spectra, they have 
performed a thorough starburst-AGN decomposition for 
each object which allowed them to estimate the fractional 
AGN contribution to the total IR luminosity output of 
each source. They demonstrate that the 4 galaxies hav- 
ing the largest mid-IR AGN fractions (> 60%) in their 
sample have q values consistent with the canonical value. 
Furthermore, they find that the FIR/radio ratio shows 
no trend with the fractional contribution of AGN activity 
in the galaxies in the IR, consistent with our results. 

6. SUMMARY AND CONCLUSIONS 

Based on a new radio-optical-IR catalog we have sep- 
arated our radio (NVSS) and IR (IRAS) detected SDSS 
galaxies (0.04 < z < 0.2) into star forming, composite, 
Seyfert, LINER, absorption line galaxies and quasars, 
and we have performed an in-depth analysis of the radio- 
FIR correlation for various types of star forming and 
AGN-bearing galaxies. Utilizing our NUV-NIR SED 
based star formation rates in combination with FIR and 
radio luminosity (expected to directly trace star forma- 
tion), we have statistically quantified the source of radio 
and FIR emission in the galaxies in our sample. We find 
that Seyfert galaxies and quasars have FIR/radio ratios 
lower than the canonical value for star forming galax- 
ies. This is due to an additional contribution to their 
radio continuum emission, which likely arises from their 
AGN. We further show that FIR-detected absorption 
and LINER galaxies are on average strongly dominated 
by AGN activity in both their FIR and radio emission; 
however their average FIR/radio ratio is consistent with 
that expected for star forming galaxies. In summary, 
our results imply that most AGN-containing galaxies in 
our sample have FIR/radio flux ratios indistinguishable 
from those of the star-forming galaxies. Thus, attempts 
to separate AGNs from star-forming galaxies by their 
FIR/radio flux ratios alone is a poor discriminant of AGN 
activity, except for the most powerful radio-loud AGN. 
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TABLE 1 

Sample summary 





IRAS (FSC + PSC) 


SDSS (MAIN + QUASAR) 


IRAS - SDSS 


total radio sample 


18313 


9591 


524 


Quasars 




4490 


21 


Absorption 




3072 


16 


Composite 




654 


203 


SF unambiguous 




621 


216 


SF ambiguous 




9 





AGN unambiguous 




454 


43 


AGN ambiguous 




291 


25 


Seyfert unambiguous 




200 


37 


LINER unambiguous 




254 


6 



The first column denotes the number of radio - IRAS (Point Source, PS, and Faint Source, FS) catalog with high quality IR photometry. 
The second column shows the number of sources in the radio - SDSS ("main" and quasar) catalog, and the third column is the matched 
radio-SDSS ("main" and quasar)-IRAS catalog. The rows indicate the various galaxy types we separate the objects into. The unambigu- 
ous/ambiguous selection is based on various spectroscopic diagnostic tools (see Fig. [5]and text for details). The shown numbers are limited 
to the 0.04 < z < 0.3 rcdshift range. 
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TABLE 2 

Fractions of star formation and AGN activity in radio and FIR regimes for various types of galaxies 





RADIO 


FIR 


SF 


AGN 


SF 


AGN 


SF 


100% 


0% 


100% 


0% 


Composite 


81.3% 


18.7% 


90.7% 


9.3% 


Seyfert 


56.8% 


43.2% 


76.1% 


23.9% 


Abso+LINER 


11.3% 


88.7% 


12.8% 


87.2% 



Fractional star formation/ AGN contribution to the radio and FIR for different optically selected galaxy types. For each population the 
fractions were estimated in a statistical man ner based on the distribution of the difference between a galaxy's SED-derived star formation 
rate and radio/FIR luminosity (see Fig. 1121 and text for details). 
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Fig. 1.— Distribution of distances between the radio and FIR detections for the NVSS-IRAS (full black line) and NVSS-SDSS-IRAS 
(dashed red line) samples. The cumulative distribution is shown in the inset. 
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Fig. 2. — The distribution of flux density at 20 cm (top panel) and 60 fim (bottom panel) for various radio-selected samples indicated in 
the top right of the panels. 
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SDSS-NVSS 




0.0 0.1 0.2 0.3 0.4 



z 

Fig. 3. — Top panel shows the 1.4GHz luminosity as a function of redshift for the NVSS - SDSS galaxies. The bottom panel shows their 
absolute optical r band magnitude (not K-corrected) as a function of redshift. 
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SDSS-NVSS-IRAS 
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z 

Fig. 4. — The top two panels are the same as Fig.[3]but for NVSS - SDSS - IRAS galaxies. The bottom panel shows the FIR luminosity 
vs. redshift. 
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log([Nll]/Ha) log([Sll]/Ha) log([Ol]/Ha) 



Fig. 5. — Optical spectroscopic diagnostic diagrams (sec Kauffmarm et al. 2003b; Kcwlcy ct al. 2001, 2006) that separate emission-line 
galaxies into star forming, composite galaxies, and various types of AGN (Seyferts and LINERs). The top panel shows the SDSS-NVSS 
sample, and the bottom panel the SDSS-NVSS-IRAS galaxies. Large symbols represent unambiguously identified galaxies (see text for 
details). Blue filled squares represent SF galaxies and green dots show composites. Red open squares, and circles represent unambiguous 
Seyferts and LINER-s, respectively. 
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Fig. 6. — Redshift distribution of galaxies in the NVSS-SDSS (top panel) and NVSS-SDSS-IRAS (middle panel) samples. The bottom 
panel shows the IR detection fraction as a function of redshift. Various galaxy types are indicated in the top right of the middle panel. 
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Fig. 7. — The radio-FIR, correlation shown for the NVSS-SDSS-IRAS galaxies (dots) and quasars (open circles). The middle panels show 
the FIR/radio ratio (i.e. q parameter defined by eq. [TJ as a function of FIR (left) and radio (right) luminosities. Note a decrease in the 
q- value with increasing radio luminosity. The bottom panels show q as a function of redshift (left) and the distribution of q with its best 
fit Gaussian (right). The number of objects in the sample (N), the biweighted mean FIR/radio ratio (< q >) and the dispersion (a) are 
indicated in the bottom right panel. 
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0.04 < z < 0.30 





Fig. 8. — Equivalent to Fig. [7] but for different galaxy types. The different symbols are indicated in the top right panel, and lines of 
constant q have been added to the top left panel. 
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Fig. 9. — The radio-FIR correlation shown for the NVSS-SDSS-IRAS galaxies divided (from top to bottom) into a) absorption line 
galaxies, b) LINERs, c) Seyferts, d) composite, and e) star forming galaxies. The left panels show the FIR-radio correlation slope, q, as 
a function of redshift. The right panels show the distribution of q for each class of galaxies in the redshift range 0.04 < z < 0.3, free of 
selection effects due to SDSS fiber-sizes used for their optical spectroscopy. 
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Fig. 10. — FIR/radio ratio (q) as a function of redshift (left panel) and the distribution of q for 21 IR (IRAS) detected quasars (right 
panel) in our radio-IR-optical sample. Note that the median q value for qauasars is lower than the average q obtained for star forming 
galaxies, but comparable to that obtained for Seyfert galaxies (see Fig. |9}. The source with q ~ —1 is a strong radio source. 
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Fig. 11. — 1.4 GHz radio (left and middle panel) and FIR (right panel) luminosity as a function of star formation rate, derived via SED 
fitting to the NUV-NIR SDSS photometry (see text for details). In the three large panels different galaxy types (symbols) and samples 
(indicated above and in each panel) are shown. The grey-scale histogram in each plot shows the distribution of star forming galaxies for 
a given sample. Superimposed on the plots (dashed lines) are calibrations converting radio and FIR luminosity to star formation rate 
HYun et 3.1.1120011: IKennicutFIITQQl) . 
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Fig. 12. — FIR excess (AlogLpm) vs. radio excess (A log Li.4GHz) for different types of NVSS-SDSS-IRAS galaxies (labeled in each 
panel) within a redshift range of 0.04 to 0.3. The FIR and radio excess has been defined a s th e difference in logarithm of the observed 
emission and that expected from their SFR, as given by the SED modeling (see eq. [5]and Fig. 1111 for details). Histograms on the axes show 
the distributions of the intrinsic AlogL values. A Gaussian fit to the distribution of star forming galaxies is also shown in each panel (solid 
line). In each panel filled dots represent the galaxies, and the gray-scale shows the distribution of star forming galaxies to guide the eye. 
Lines of constant q are also shown. 
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